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Abstract. The present status of lattice calculations of /b, Bb and SU(3) flavour breaking ratios such 
as /b s //b is reviewed. Particular attention is devoted to systematic uncertainties, such as those arising 
from the lack of simulation data for dynamical quarks with realistic masses, and the related difficulties 
associated with chiral extrapolations. Global averages for decay constants and mixing parameters are 
presented (Tabled, and the procedures to obtain them are discussed in detail. 
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1 Introduction 

A lot of activity is currently devoted to pin down the ele- 
ments of the CKM matrix Vckm- ln the Standard Model 
Vckm is unitary, which implies triangle relations like 



X 



VudV* b + v cd v; h + v td v t i = o. 



(i) 



Any deviation from unitarity is interpreted as a signa- 
ture of "new physics". In order to probe this scenario, ex- 
perimental and theoretical inputs are being used to over- 
constrain the elements of Vckm- However, relations be- 
tween measurable quantities and CKM matrix elements 
involving heavy quarks are usually afflicted with large 
hadronic uncertainties. A typical example are the mass 
differences AM d and AM S in the B° B° and B° B° sys- 
tems: 
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The limited accuracy with which the CKM elements on 
the rhs. are known comes from theoretical uncertainties in 
the decay constants /b, /b s and the B-parameters -Bb and 
Bb s - These quantities have been computed using lattice 
calculations, an approach which was specifically designed 
for a systematic non-perturbative treatment of QCD. In 
order to chart the progress made and to provide global 
estimates for these quantities, the CKM-Lattice Working 
Group was founded in February 2002 PE]. Here I report 
on recent results and present global averages. 
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2 Heavy quarks on the lattice 



Systematic effects in lattice simulations, as well as the 
more specific problems of treating heavy quarks on the 
lattice have been described many times in the literature 
(see e.g. PITT]). 

The great majority of lattice results for heavy-light 
decay constants and B-parameters have to date been ob- 
tained in the quenched approximation, where quark loops 
are neglected in the evaluation of observables. The effects 
of non-zero lattice spacing a (lattice artefacts) have been 
studied extensively, though, and in many cases an extrap- 
olation to the continuum limit was performed. In order 
to guarantee a smooth continuum behaviour, the non- 
perturbative renormalisation of quark bilinears and four- 
fermion operators in the discretised theory proved to be 
instrumental 5,6]. For some quantities, the level of pre- 
cision that can be reached in the continuum limit in the 
quenched approximation is about 5%. This then implies 
that current estimates for decay constants are completely 
dominated by quenching effects. As a consequence, most 
collaborations now focus on simulations with Nf = 2 or 3 
dynamical quark flavours. 

Another important systematic effect is the use of un- 
physical values of the light quark masses m u , m d , both in 
quenched and unquenched simulations. In particular, com- 
monly used algorithms for dynamical quarks slow down 
considerably for masses smaller than about m s /2. Several 
proposals to address this problem have been made 0ODIII1 
ITT)] , but the effectiveness of each approach must be studied 
in more detail. In order to make contact with the chiral 
regime, one currently relies on extrapolations in the light 
quark mass, using Chiral Perturbation Theory (ChPT) as 
a guide. It has only been realised relatively recently that 
this can introduce large uncertainties, since chiral loga- 
rithms are not necessarily under control jl lU12lll.{j . 
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Whenever one deals with heavy quarks on a lattice of 
spatial extent L and lattice spacing a, one is faced with a 
multi-scale problem, in the sense that the following three 
inequalities cannot be satisfied simultaneously: 

am h < 1, m^L > 1, L/a < 50. (4) 

Violation of the first relation implies the presence of large 
lattice artefacts, the second inequality must be satisfied if 
one aims for small finite-volume effects, and the third is 
dictated by capacities of current computers. 

Several strategies to deal with this problem have been 
applied over many years, among them the "static approx- 
imation" [T"I], the non-relativistic formulation (NRQCD) 
"T"""| , the so-called "Fermilab-approach" ,]*""J , and extrapo- 
lations in the heavy quark mass. 



3 Recent results 

We now discuss results from several recent calculations, 
which will also illustrate some of the issues raised earlier. 

The first topic is a recent benchmark calculation of 
/d s in the quenched approximation |17| . For lattice simu- 
lations, the D s meson is particularly appealing, since both 
the charm and the strange quark can be treated directly, 
i.e. no extrapolations are required to make contact with 
the physical values of the valence quark masses. In ref. 
"T7) the potentially large lattice artefacts arising from rel- 
ativistic c-quarks are eliminated through an extrapolation 
to the continuum limit. By means of employing a lattice 
action in which the leading lattice artefacts of 0(a) were 
removed non-perturbatively [IB|, the convergence of the 
results to the continuum could be accelerated. Further- 
more, the authors used a non-perturbative estimate of the 
renormalisation factor which connects the axial current on 
the lattice with its continuum counterpart |1 9| . 

The main result in ref. |17| is the continuum result 
for /d s in units of the hadronic radius ro 20 , namely 
'"o/Ds = 0.638 ± 0.024, which translates into 

/ Ds = 252 ± 9 MeV (5) 

if the phenomenological value ro = 0.5 fm is inserted. It is 
worth emphasising that the only remaining uncertainty is 
due to quenching. A crude estimate of the quenching error 
is obtained through the scale uncertainty, i.e. the fact that 
different quantities yield different estimates for the lattice 
scale in the quenched approximation. The typical size of 
this ambiguity is about 10% |21| . 

In another recent calculation by de Divitiis et aL[22] a 
new strategy to deal with the multi-scale problem was ap- 
plied. Here the condition m^L 1 was sacrificed in favour 
of amb <C 1. In this way one is able to accommodate a fully 
relativistic 6-quark, and a physically meaningful result is 
obtained if one succeeds in determining the distortion due 
to the unphysically small volume. The key observation is 
that this can be achieved through a sequence of finite-size 
scaling steps, which relate the results obtained for several 



lattice sizes Lq, L\ 1 . . .: 

MLo) -> MLi) -> M**) -> ••■ (6) 
f B {L k+ i) = fs{L k ) cr B (L k ), L k+ i > L k (7) 

Since the value of /b for the smallest volume, i.e. /b(-^o) 
and the so-called "step-scaling function" a-Q{L k ) can be 
computed in the continuum limit, one obtains a result for 
/b for physically large volumes in a controlled manner. 
The main assumption, which was verified explicitly in ref. 
|22j . is that the finite-size effects depend only weakly on 
the heavy quark mass. Starting from Lq = 0.4 fm, a total 
of two scaling steps were performed, each time doubling 
the lattice size. The initial value of the decay constant in 
the continuum is Jb(Lo) = 471(2) MeV, and together with 
the continuum estimates for the step scaling functions, 
cr B (L ) = 0.400(3) and ob(Li) = 0.92(4) one obtains 

/ B = 173±8±4MeV, (8) 

at L2 — 1.6 fm, which is large enough to be identified with 
the infinite volume limit. The quoted uncertainties cover 
all errors, except that due to quenching. In a similar way 
the authors of ref. [221 obtain 

/ Ba = 194±6±4MeV, (9) 
/ D = 217±7±5MeV, / Ds = 239 ± 5 ± 5 MeV, 

the result for /rj s being compatible with ref. |17j . 

As mentioned earlier, the issue of chiral logarithms 
in SU(3)-flavour-breaking ratios such as /b h //b or £ has 
attracted a lot of attention recently. The dependence of 
heavy-light decay constants on the mass of the light quark 
has usually been modelled according to a naive linear 
ansatz, in which chiral logarithms, as well as analytic 
terms arising at NLO in ChPT, are neglected. The full 
expression at NLO reads 

-Bl _ 1 = ( m 2 _ m 2)J 2 ( M ) 
IB 

[iM^k) + |/pK) - fipKO] (io) 

where ip(mps) = rap s ln(mp s //i 2 ) and /2 is a low-energy 
constant. As was pointed out by Kronfeld & Ryan ["Q], the 
inclusion of chiral logarithms in the chiral extrapolation of 
lattice data for heavy-light decay constants can drastically 
change /b s //b and consequently also £, which enters fits 
to the CKM parameters. By assuming g 2 — g 2 D * D7T — 0.35 
[2"3"| and fi = 0.5(3) GeV -2 , Kronfeld and Ryan conclude 
that £ = 1.32 ± 0.10, which is more than 10% larger than 
the global estimate quoted by Ryan in 2001 [2I|. It is some- 
what ironic that the quantity £, which for a long time had 
been assumed to be only weakly sensitive to systematic ef- 
fects, should be subjected to such a large uncertainty. Note 
that the corresponding ratio of B-parameters, Bb s /Bb is 
largely unaffected, since the coefficient of the chiral loga- 
rithm is oc (1 — 3<7 2 ), which is close to zero. 

The predicted enhancement of /b s //b due to chiral 
logs seems plausible, since the corresponding ratio in the 



3ii decays and mixing 



3 




0.0 2.0 4.0 6.0 

( r <rJ 2 

Fig. 1. Dependence of <Pf B and <&f Bis on the light quark mass. 

light quark sector, fn/ fn, is known to come out too small 
in quenched QCD, if only a naive linear quark mass depen- 
dence is assumed [25] . In ref. ^2! it was argued that the 
chiral logarithms in /b b //b and /k//tt are nearly of the 
same size, so that one would expect /b b //b ~ /k/ fn = 
1.22, which is compatible with 111) , but closer to previous 
global estimates. The key question for any future deter- 
mination is whether or not the quark masses used in the 
simulation are light enough so that ChPT at NLO gives a 
good description. In the context of /b b //b this was stud- 
ied recently by the JLQCD Collaboration 26 . The depen- 
dence of <£/ B = /b\/A/b and <Pf Bs on the light quark mass 
is shown in Fig.^ JLQCD conclude that chiral logarithms 
are not observed in the studied mass range. The modelling 
of the effect due to chiral logs yields a drop of up to 11% 
in /b relative to /b b . The corresponding enhancement of 
/b / /b is of the same order of magnitude as that of ref. 

ED 



4 Global estimates: Strategies & Results 

Before I present global averages for decay constants and 
£?-paramctcrs, I would like to outline the strategy which 
was used to obtain them. The main point to note is that 
a global analysis of lattice results is complicated by the 
fact that systematics can vary substantially among results 
from different collaborations. The main differences lie in 
the treatment of lattice artefacts, the choice of quantity 
that sets the lattice scale, the renormalisation procedure 
for currents and four-fermion operators on the lattice, and 
the details of extrapolations in the quark masses, either 
down to the chiral regime, or, where applicable, to the 
mass of the &-quark. Therefore, a straightforward global 
average of results is unreliable. 

For the procedure applied here I have decided to focus 
on the quenched approximation: although unphysical, the 
quenched approximation is known to describe the light 
hadron spectrum at the level of 10% [23- Furthermore, a 
wealth of results is available for Nf — 0: the continuum 
limit has been taken in almost all recent quenched cal- 



culations, and other systematic effects have been studied 
thoroughly. 

Simulations with dynamical quarks, either for Nf = 2 
or 3 have not yet reached the same level of maturity. The 
masses of the sea quarks are still quite large, so that their 
effects on observables are likely to be suppressed. Perform- 
ing reliable continuum extrapolations is much more ex- 
pensive, all of which implies that a clear separation of sea 
quark effects from lattice artefacts is very difficult. There 
are, however, attempts to expose the effects of unquench- 
ing by focusing on ratios in which systematic uncertain- 
ties other than quenching largely cancel. For instance, the 
CP-PACS E3 and MILC M collaborations have found 

fSr 2 /f^r° ~ L1 ° at a ~ °- lfm > which implies a 10%- 
enhancement in /d b due to dynamical quark effects. No 
extrapolations in the valence quark masses are required, 
but it remains to be seen whether this number is stable 
against lattice artefacts and whether lighter sea quarks 
lead to a significant change. 

As regards chiral logarithms in ratios like /b s //b and 
£, their effects have not been quantified from first princi- 
ples so far. In what follows, I shall therefore use results 
assuming the naive (LO) dependence on the light quark 
mass, but allow for a 10% increase in /b.//b and /d b //d- 

The starting point for the global analysis is /d b in the 
quenched approximation ( Nf — 0) . The recent benchmark 
calculation [T7 demonstrates that the continuum value in 
quenched QCD (for a given quantity that sets the scale) 
is obtained with a total accuracy of 4%. The central value 
may vary between 225 and 255 MeV, depending on the 
chosen scale. Thus we have 

/^ t=0 = 240 ± 10 (stat) ± 15 (scale) MeV. (11) 

In order to estimate /d one can divide this result by 
f D Jf D = 1.12 ± 0.02^°;^. This is the number quoted 
by Ryan for Nf — 0, except for the additional 10% 
asymmetric uncertainty due to chiral logs. This yields 

/j^=° = 214 ± 10 (stat) ± 13 (scale)t^ (xlog) MeV. (12) 

In Fig.|2]the quenched global estimates of eqs. (|llf) and (|12f> 

are compared with recent determinations from various 
groups, who have used either the Fermilab approach |29l 
l3Till31ll27ll28| . or the O(a) improved Wilson action [3Hl3"3"l 
I34H17U22| to treat the heavy quark. Note that the results 
from different collaborations have not been converted to 
a common lattice scale. Without this conversion it is not 
clear whether the observed spread is just indicative of the 
scale ambiguity, or due to some other, possibly uncon- 
trolled systematic effect. 

The next step is to account for dynamical quarks. In 
my view the safest procedure is to multiply the quenched 
values of /d and /d b by the ratio 

fp^ =2 /fp^ = ° = 1.10±0.05, P = D,B, q = d, s, (13) 

in which some of the systematic errors can be expected to 
cancel. The number in eq. I|13l) is based on observations 
in ref. |27| that dynamical quarks enhance the values of 
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Fig. 2. Recent results for decay constants and B-parameters in the quenched approximation, compared to the global averages 
of eqs. ITU . fl3 ) CH , and lO (solid vertical lines). For D -mesons the dashed lines indicate the statistical error, while 
the systematic error is represented by the dotted lines. For the B-system the dashed vertical lines denote the total error, with 
the uncertainty due to chiral logs included in quadrature. Dotted error bars on individual determinations indicate the quoted 
systematic error. The result for Bb and £ in |2G| was obtained for Nf = 2. 



decay constants by 5 — 15 %, largely independent of the 
valence quark contents. Thus we find 

fg*= 2 = 264 ± 11 (stat) ± 22 (quen) MeV (14) 
/^f= 2 = 235 ± 11 (stat) ± 19 (quen) t°i ( X log) MeV, 

where I have combined the scale uncertainty and the error 
in eq. Ijl3(l into a total quenching error. Here we make no 
attempt to try and estimate the effects of a dynamical 
strange quark, and hence the above numbers are our final 
results for /d and /d s which we list once more in Tabled 
The procedure to obtain estimates for /b and /b s is 
entirely analogous, with one important difference: unlike 
the case of D-mesons, there is not yet a quenched bench- 
mark calculation of /b s , which would yield the total error 
within the quenched approximation. As a starting point 
we therefore use the "global representation" of quenched 
results quoted by Ryan [23], i.e. 

/^ f=0 = 200 ±20 MeV. (15) 

Here the error has been obtained by requiring consistency 
with a number of different results being subjected to dif- 
ferent systematics. Dividing by [23] 

/B 8 // B = 1.15±0.03i^( x log) (16) 

gives 

/^=° = 174 ± 18 ± x ° ( X log) MeV, (17) 

and multiplication with the ratio of eq. (|13fl yields the 
numbers listed in Tabled In Fig.[2]the global quenched es- 
timates are again compared to individual results obtained 
using NRQCD |331l3^1l37| , the Fermilab approach [2H5U1 
I31U27U28) . the finite-size scaling technique [22] , as well as 
extrapolations in the heavy quark mass from the region of 
the charm quark mass | 38U33U34 | . 



Table 1. Global estimates for decay constants and B- 
parameters. The numbers are understood to refer to full QCD. 

/ Da = 264 ± 11 ± 22 (quen) MeV 

/d = 235 ± 11 ± 19 (quen) ± 2 ° (xlog) MeV 

fa, = 220 ± 25 MeV 

/b =191 ±231^ ( X log) MeV 

Bb = Bb s = 1.34 ± 0.12 

£ = 1.15 ± 0.05 ±°;Jg 

/b b \/Bb7 = 255 ± 31 MeV 

/b \/Bb~ = 221 ± 28 1 2 °2 (xlog) MeV 



For the -B-parameters Bb, Bb s and the ratio £, there 
are not so many results available, and only one recent cal- 
culation uses dynamical quarks |26j . It is then clear that 
systematics cannot be studied as thoroughly as for decay 
constants. In particular, the continuum limit has not been 
taken in any study so far. It turns out, though, that results 
for £>-parameters and £ are broadly consistent, regardless 
of whether NRQCD [39,40,26 , relativistic heavy quarks 
34,38,41 or the static approximation 42 43, 44, 45 are 
used in simulations. Apparently, dynamical quarks do not 
lead to a significant enhancement of Bb or Bb s |2H] > con- 
trary to what is observed for decay constants. Chiral log- 
arithms in the SU(3)-flavour breaking ratio Bb s /Bb are 
suppressed, and a linear extrapolation in the light quark 
mass yields a value compatible with one. All published 
data are then consistent with the global representation: 

B B (m b ) =0.85 ±0.08 => B^ LO = 1.34 ± 0.12 
B Ba (m b ) = 0.85 ±0.08 =*> B^° = 1.34 ± 0.12 (18) 
BbJBb = 1.00 ± 0.03, £ = 1.15 ± 0.05±°;^ ( x l og ). 
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The value of £ is obtained by combining eq. (|16Jl with 
the above result for BbJB b . The collection of global es- 
timates in Tabled is consistent with refs. [Tll24ll47| . 

5 Future studies 

How can the current global estimates be improved in order 
to sharpen the constraints on CKM parameters? There are 
several areas in which progress can be made. 

Future efforts must surely focus on dynamical simu- 
lations, but the quenched approximation remains help- 
ful for the understanding of several issues. In my view, 
a quenched benchmark value of /b b is very important to 
quantify the uncertainties arising from discretisation and 
renormalisation effects. Although this is precisely the aim 
of the finite-size scaling method of [22], an independent 
check using a different method should be performed. Here 
one may think of an interpolation between results ob- 
tained in the static approximation and those obtained near 
m c , after the continuum limit has been taken. With the 
advent of methods that allow for precise non-perturbative 
determinations of the renormalisation factor for the axial 
current [JB], as well as for a better signal/noise ratio in 
the static approximation |49| . one can obtain results with 
much higher accuracy than previously possible jSU]. This 
strategy may also be extended to B-parameters |51) . 

Obviously the most pressing problem is to simulate 
small dynamical quark masses more efficiently. The role of 
"improved staggered quarks" has been vigorously empha- 
sised [5] in this context, with first results being reviewed 
in [S^J. Eventually such efforts should yield not only es- 
timates for, say, /^ f_3 , but should also settle the issue 
of chiral logarithms in /b,//b, £ and /k//tt- The latter 
is, in fact, explored separately by a number of groups us- 
ing different algorithms and discretisations |5.'iU26| . In this 
context, the so-called Grinstein ratio (/b b //b) /(/d s //d) 
in which the chiral logarithms cancel [53], has been sug- 
gested as a more reliable way to determine /b b / /b on the 
lattice |2E], once /d b and /d have been measured at high- 
luminosity charm factories such as CESR-c [55]. This ap- 
proach, however, can only succeed if the quark masses 
used in simulations are light enough so that the data are 
consistently described by ChPT at NLO. 

Finally, one may think of better ways to obtain global 
averages. This may entail potentially laborious re-analyses 
of existing simulation data. In order to facilitate a com- 
parison all results should be converted to a common scale, 
such as ro. Systematic errors should best be quantified 
after an extrapolation to the continuum limit. Here on 
might concentrate on ratios like /b b //d b for fixed Nf, or 
/ D f_2 // D f=0 . Thereby one has a better chance to expose 
the effects of replacing c by 6-quarks and those due to 
dynamical quarks. 

It is a pleasure to thank P. Hernandez, M. Papinutto 
and C. Pena for helpful discussions. I thank the organisers 
as well as the conveners V. Lubicz and B. Pioline for a 
stimulating conference. I am grateful to I. Montvay and 
C. Pena for a critical reading of the manuscript. 
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